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Abstract 


An integrated thermal -structural finite 
element approach for efficient coupling of 
transient thermal and structural analysis is 
presented. New integrated thermal -structural rod 
and one dimensional axisymmetric elements 
considering conduction and convection are 
developed and used in transient thermal -structural 
applications. The Improved accuracy of the 
integrated approach is illustrated by comparisons 
with exact transient heat conduction-elasticity 
solutions and conventional finite element 
thermal -finite element structural analyses. 

Results Indicate that the approach offers 
significant potential for further development with 
other elements. 
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Nomenclature 

Inner radius of cylinder, see 
fig. 4 

cross-sectional area 
outer radius of cylinder, see 
Hg. 4 

st ra 1 n -d 1 s p 1 acement 
Interpolation matrix 
temperature gradient 
interpolation matrix 
specific heat 

arbitrary constants, see equation 

(8) 

finite element capacitance matrix 
elasticity matrix 

forcing function, see equation (7) 
homogeneous solution, see equation 

. 

particular solution, see equation 

( 8 ) 

modulus of elasticity 
force 

finite element nodal force vector 
convective heat transfer 
coefficient 
thermal conductivity 
thermal conductivity matrix 
finite element stiffness matrix 
finite element conductance matrix 
for conduction 

finite element conductance matrix 
for convection 
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Subscripts 

c 

e 

h 

S 

T 

Superscript 

T 


length 

conduction-con vection rod 

parameter, m ^hp/kA' 

order of diff^entlal equation, 

see equation (7) 

finite element Interpolation 

function for nodeless variables, 

see equation (10) 

finite element Interpolation 

functions 

finite element displacement 
Interpolation functions 
finite element temperature 
interpolation functions 
perimeter 

surface heating rates 
radial coordinate, see fig. 4 
volumetric heat generation rate 
finite element heat load vector 
finite element surface area 
time 

temperature 
nodeless parameter 
nodeless variable 
reference temperature for zero 
stress 

environmental temperature for 
convective heat exchange 
displacement components 
ln(b/a), see table 4 
finite element volume 
cartesian coordinates 
nondlmensional coordinate, X = x/L 
vector of thermal expansion 
coefficients 
density 

dependent variable In differential 
equation, see equation (7) 


conduction heat transfer 
element matrix or vector 
convective heat transfer 
structural 
thermal 


transpose of a matrix 
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A recent overview^ of research In structures 
and materials for future space transportation 
systems cited the challenge of providing 
lightweight structures for large space 
transportation vehicles that are repeatedly 
subjected to severe aerodynamic heating and yet 
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must have structural efficiency, reliability and 
durability approaching that of commercial 
aircraft* One of the most challenging 
thermal -structural analyses faced by aerospace 
engineers today is the Space Shuttle Orbiter. 
Shuttle design experience has shown that the 
design of large structures for hostile reentry 
environment taxes the capability of existing 
thermal and structural analysis methods. The 
complex thermal behavior of the Shuttle has 
required a large number of highly detailed lumped 
parameter analytical models. A review^ of needs 
revealed by the Shuttle thermal analysis cited the 
need for improved modeling procedures to reduce 
model size and more efficient methods of 
transferring data between thermal and structural 
models. Experience in the repetitious thermal and 
structural analyses required for the design of 
convectively-cooled structures and in the optimum 
design^ of structures subjected to heating and 
mechanical loads at elevated temperatures also 
points to the need for more efficient coupling of 
the thermal and structural analyses. 

Thermal analysis of complex structures is 
generally carried out by one of two alternative 
methods: (1) the finite difference lumped 

parameter method, or (2) the finite element 
method. The lumped parameter approach is most 
widely used, but the finite element method is an 
attractive alternative since it provides 
capabilities for both thermal and structural 
analysis of general structures. Frequently the 
thermal analysis is performed by the lumped 
parameter method and the structural analysis by 
the finite element method. Because of basic 
differences between the analytical models an 
efficient interface is difficult to achieve. 

The finite element method offers the greatest 
potential for efficient coupling of the thermal 
and structural analyses, but the historically 
better capabilities and efficiency of the lumped 
parameter method compared to early finite element 
methods has so far prevented the finite element 
method from receiving widespread acceptance for 
thermal analysis. Recent experience^ with the 
finite element method for combined conduction 
forced convection analysis has shown that: (1) 

the methods currently have about the same analysis 
capabilities, (2) the finite element method has 
superior accuracy, and (3) the lumped parameter 
method maintains an edge in efficiency, but that 
as additional finite element experience is gained, 
improvements in efficiency can be achieved. 

Research programs are currently underway at 
Langley Research Center to improve both the 
capabilities and efficiency of the finite element 
thermal analysis method and develop more efficient 
coupling between the finite element thermal and 
structural analyses. The authors previously 
presented a paper^ which focused on finite element 
methodology for efficient coupling of steady-state 
thermal and structural analyses. The purpose of 
this paper is to extend the methodology, denoted 
as integrated thermal -structural analysis, for the 
transient -thermal static-structural analysis case. 

Characteristics of Integrated 
thermal -structural analysis are first discussed. 

A nodeless variable approach for the development 
of transient integrated thermal -structural 
elements is then described. Next, the nodeless 


variable approach is used to develop new 
Integrated thermal -structural rod and 
one-dimensional axisymmetric elements. Finally, 
the accuracy and efficiency of the integrated 
approach is demonstrated by solving three 
thermal -structural examples by both the 
conventional and integrated finite element 
approaches. In two of the examples, the accuracy 
of the finite element solutions is evaluated by 
comparisons with exact transient heat 
conduction-elasticity solutions. Finite element 
matrices are presented in Appendices. 


Integrated Thermal -Structural Analysis 
Finite Element Analysis 

Finite element (F.E.) formulations for 
transient thermal problems are derived by the 
method of weighted residuals. ° In general, 
element temperature T(x,y,z,t) and temperature 
gradients are expressed in the form 

{T}= [NT]{T(t))e (]a) 

(S/ayl = CBT]{T}e (lb) 

(aT/sz) 

where [Nj] denotes a matrix of the temperature 
interpolation functions, [Bj] denotes a matrix 
of temperature gradient interpolation functions 
and {T(t)}g denotes a vector (one dimensional 
array) of nodal temperatures. For transient 
thermal analysis of conduction with convection 
boundary conditions the equations for a typical 


element are 


[cDgftJe + [Kc + 

(2) 

where [cL is the element capacitance matrix, 

[Kc^e 3nd [K^L are element conduction and 
convection conductance matrices, respecti vely, and 
{Q}g is a vector of nodal heat loads. These 
matrices. In general, are expressed in the form of 
integrals over the volume, Vg, and surface, Sg» of 
an element. The element equations are 

[CL =/ PcENjl^CNTldV 

(3a) 

c 

[KcL =/ [Bj3'''Ck][BT]dV 

(3b) 

w 

[KhL =/ b[Nx]T[NT]dS 

(3c) 

■CQ^e Q[Nj]^dV + / q[Nj]^dS 

^e ^e 

(3d) 


+ / h T„CNj]''’dS 
^e 

where p is the density, c is the specific heat, 

[k] denotes the conductivity matrix, h is a 
convection coefficient, Q is an internal heat 
generation rate per unit volume, q is a surface 
heating rate per unit area and T„, denotes the 
environmental temperature for the convective heat 
exchange. The superscript T denotes the transpose 
of a matrix. The convection coefficient, the 
internal heat generation rate, the surface heating 
rate and the environmental temperature are, in 


2 



general, time-dependent* All thermal parameters, 
herein, are assumed constant although they may be 
temperature dependent in general. Finite element 
formulations for other heat transfer modes such as 
forced convection and radiation are given in 
references 4 and 7, respectively. The simultaneous 
differential equations (2) are typically solved by 
time-marching schemes so that the temperatures {T} 
are computed at discrete time values. 

In general formulations of transient 
thermal -stress problems, the heat transfer and 
elasticity solutions are coupled through a 
mechanical coupling term in the heat conduction 
equation and inertia terms in the elasticity 
equations. In most engineering applications the 
coupling effects are not significant, and the 
coupling terms are neglected. In the uncoupled 
formulation, transient temperatures are computed in 
the thermal analysis and are used as input to a 
quasi-static structural analysis. Finite element 
formulations for the structural analysis are 
usually derived from a variational formulation.® 
Since the structural temperatures vary with time, 
the structural analysis consists of a sequence of 
static analyses at selected time values in the 
transient thermal analysis. In a typical 
structural analysis, element displacements (u,v,w) 
are expressed in the form 

{u} = [Ns3{u}e (4) 

where [N5(x,y,z)] denotes the displacement 
interpolation functions, and {u}g denotes a vector 
of element nodal displacements. For the structural 
analysis the equations for a typical element are 

[K]e{u}g - {F>e (5) 

where [K]g is the element stiffness matrix, and {F} 
is a vector of nodal forces which can consist of 
mechanical forces and/or equivalent thermal forces. 
Only the equations for the equivalent thermal force 
will be shown herein. The element equations are 


[K]e =/ [Bs]T [D][Bs]dV 
^e 

(6a) 

{F>e =/ [65]"'’ [D]{a}T dV 

(6b) 


e 


where [B5] is the strain-displacement 
interpolation matrix, [D] denotes the elasticity 
matrix, and {a} is a vector of thermal expansion 
coefficients. T(x,y,z,t) is the temperature within 
the element computed in the transient thermal 
analysis. 

The temperatures enter the structural analysis in 
two ways: (1) the structural elasticity matrix [D] 

and thermal expansion coefficient vector {a} are, 
in general, temperature dependent, and (2) the 
equivalent thermal forces, equation (6b), involves 
an integration of the temperature T(x,y,z,t) over 
the element volume. For a temperature dependent 
elasticity matrix, the element stiffness matrix, 
equation (6a), varies during the transient 
response. In this paper the elastic properties and 
the thermal expansion coefficient are assumed 
constant, hence the element stiffness matrix is 
constant and is computed only once. The nodal 
forcevector is computed for each temperature vector 
at the selected time values, and consequently the 
structural analysis is a static analysis with 
multiple load cases. 


Integrated Analysis 

A step toward integrated thermal -structural 
analysis capability is the use of a common 
methodology in a single program such as NASTRAN® 
which contains both thermal and structural analysis 
capability. One disadvantage of this approach, 
herein called the conventional approach, is that 
the transfer of data between analyses is often 
Inefficient because of inherent differences between 
the thermal and structural models. Another 
disadvantage of the conventional approach is that 
basic differences between the thermal and 
structural analysis requirements are not recognized 
and exploited. For example, the very common beam, 
plate and shell structural finite elements have no 
thermal counterparts in typical finite element 
thermal analysis programs because conventional 
finite elements have no provision for computing 
thickness-temperature gradients. This deficiency 
means not only that the thermal and structural 
finite element models of the same structure differ, 
but that thermal load data (e.g. temperature 
gradients) are not directly supplied by the thermal 
analysis. 

To exploit more fully the capabilities of the 
F.E. method the concept of integrated 
thermal -structural analysis was proposed in 
reference 5. An integrated thermal structural 
analysis is characterized by: (1) thermal and 

structural finite elements formulated with a common 
geometric discretization with each element 
formulated to suit the needs of their respective 
analysis, (2) thermal and structural finite 
elements which are fully compatible, and (3) 
equivalent thermal forces which are based upon the 
consistent finite element force vector computed by 
equation (6b). 

The concepts of conventional and integrated 
thermal -structural analysis approaches are compared 
schematically in Fig. 1. The sequence followed in 
a conventional analysis is shown in Fig. 1(a). 
First, a thermal analysis is performed based upon a 
thermal model selected to best represent the heat 
transfer problem. The thermal model may be based 
on the lumped parameter or F.E. method. The 
thermal analysis is followed by transfer of the 
nodal temperature data to the structural analysis. 
In most cases this data transfer involves data 
processing to bring the input temperatures in 
conformity with the needs of the structural model. 
Often the structural and thermal models use 
different nodes and elements, and approximate 
thermal forces are computed from average element 
temperatures. An integrated analysis is shown in 
Fig. 1(b). The thermal and structural analyses are 
characterized by a common model based upon thermal 
and structural elements formulated to best suit the 
respective analysis. The transfer of data is 
compatible with no data processing required. 
Consistent thermal forces are computed from thermal 
elemental and nodal input data supplied directly 
from the thermal analysis. 

Integrated Elements 

The concept of integrated transient 
thermal -structural analysis is illustrated with a 
new rod element and a new one -dimensional 
axisymmetric element based on nodeless variable 
formulations. The one-dimensional integrated 
elements consider combined conduction and 
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convection with constant thermal parameters. The 
nodeless variable approach" was utilized In 
reference 5 to develop an Integrated rod element 
for steady state analysis. The nodeless variable 
approach as used herein consists of utilizing 
Improved temperature Interpolation functions which 
employ an element or "nodeless" variable. The 
nodeless variable approach has been used 
previously In structural analysis by adding extra 
interpolation functions (sometimes called "bubble" 
modes) to the conventional interpolation 
functions. A unique feature of the nodeless 
variable approach for steady state analysis is the 
use of an exact temperature variation to compute 
the element conductance matrices, heat load 
vectors and equivalent thermal forces. The exact 
temperature variation is also used to compute an 
exact displacement variation In the structural 
element. 


Exact Steady-State Formulation 


Finite elements which yield exact values of 
the nodal variables will herein be called exact 
finite elements. Structural rod and beam elements 
under simple loadings are known to produce exact 
values of the nodal displacements. Linear 
conduction elements under simple heating also 
produce exact values of nodal temperatures." 
Further, in one-dimensional combined 
conduction-forced convection heat transfer an 
exact F.E. has been derived^ using the Galerkin 
method with upwind weighting functions. In this 
section an approach" for deriving exact finite 
elements In one-dimensional problems will be 
described. The approach Is based upon using the 
exact solution to the corresponding linear, 
steady-state differential equation to derive the 
element Interpolation functions. If the elements 
are used alone (i.e. without connections to other 
element families) an exact solution will be 
obtained. Approximate solutions will be obtained 
if the elements are used in transient or nonlinear 
applications, or if the elements are connected to 
elements of other non-exact interpolation function 
families. 


Consider a linear, nonhomogeneous 
differential equation 

ap dH(|)/dx*^ + ap^i dn-^/dx^^-l 

. X ( 7 ) 

+ ... + ao4> = f(x) 

where 0(x) is the dependent variable. The exact 
solution has the form 

't'(x) = Z) c^f.(x) + g(x) (8) 

i =1 


where c^ are arbitrary constants, f^(x) are 
typical functions in the homogeneous solution, and 
g(x) Is the particular solution. A F.E. with n 
degrees of freedom is formulated based upon the 
continuity requirements for nodal variables. 

For example, if the dependent variable only is 
required to be continuous at element nodes then an 
element with n nodes will permit an exact 
solution. Since the differential equation has n 
constants of integration the element interpolation 
functions are determined by Imposing the 
conditions 

<f'(xi) = (!>i i = 1,2 (9) 


where x^ are the nodal coordinates, and are the 
nodal values of the dependent variables. To 
accommodate the particular solution, g(x), the 
element Interpolation function is written in the 
form 

<p{x) = Njj(x)iJ)q + E N^(x),t,j (10) 

where (|)q Is a nodeless variable.® In applications 
4)q is selected as a convenient parameter and 
consequently is a known quantity a priori. Note 
that since N^-(x^*) = 1 to satisfy eq. (9), then 
Nq(x^) = 0. With the exact Interpolation 
functions known, the element matrices can be 
derived using typical finite element integrals 
such as eqs. (3) for thermal analysis or eqs. (6) 
for structural analysis. Element matrices are 
generally of size n + 1 due to the presence of the 
nodeless variable. In many cases (see Appendices 
A and B) the extra equation is uncoupled so that 
In the computer analysis only matrices of size n 
are employed. 

Transient Formulation 

In Ref. 5, the nodeless approach was used to 
develop an exact Integrated thermal -structural rod 
element for steady-state conduction and convection 
analysis. One dimensional elements such as the 
rod element are characterized by two nodal 
unknowns and one nodeless variable for the 
steady-state thermal or structural element 
formulations where the nodeless variable is 
identified as a convenient physical parameter and 
Is known for each element. 

For general transient analysis it is not 
possible to formulate closed-form Interpolation 
functions which are exact solutions since general 
transient solutions to the heat 
conduction-elasticity equations are infinite 
series. However, the nodeless variable 
interpolation function concept introduced in eq. 
(10) can be extended for the transient case to 
give accuracy superior to the conventional finite 
element approach. To extend the nodeless variable 
approach to the transient case, two approaches are 
investigated: (1) the nodeless variable approach 

previously described for steady-state analysis, 
which hereafter Is referred to as the nodeless 
parameter formulation, and (2) a nodeless variable 
approach where the nodeless variable is an unknown 
function of time. 

Nodeless parameter approach . In the nodeless 
parameter approach the temperature interpolation, 
eq. (10), for an element with two nodes has the 
form 

T(x,t) = Nq(x)To + Ni(x) Ti(t) 

+ N2(x) T2(t) 

where Tq is a known nodeless parameter, and T^(t), 
T2(t) are unknown time-dependent nodal 
temperatures. The element equations, eq. (2), for 
a typical element based on this approach may be 
written as 
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where the element subscript e has been omitted* 

The salient characteristic of this element 
formulation is that the first equation Involving 
the nodeless parameter is uncoupled from the nodal 
unknowns in the second and third equations. Thus 
the element matrices have two unknowns as for a 
conventional element. 


Nodeless variable approach. In the nodeless 
variable approach the temperature Interpolation, 
eq. (10), for an element with two nodes has the 
form 


T(x,t) =No(x) To(t) + Ni(x) Ti(t) 
+ N2 (x) T2(t) 


(13) 


where Tg(t) Is an unknown time-dependent nodeless 
variable, and Tj(t), T 2 (t) are unknown 
time-dependent nodal temperatures. Note that for 
steady-state solutions the nodeless variable has a 
value equal to Tq. The element equations, eq. 

(2), for a typical element now have the form 


CqO CqI Cq2 
ClO ^11 Ci2 
C20 C21 C22 


To 

Tl 

T2 
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ho 
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(14) 


Since the nodeless variable is unknown, the 
equations are coupled through the capacitance 
matrix due to the presence of Tq. Thus the 
element matrices have three unknowns, one more 
degree of freedom than a conventional element. 


A unique feature of the nodeless parameter 
and nodeless variable interpolation equations is 
that in the steady-state each reduce to the exact 
solution. This means the transient response can 
start from an exact initial temperature 
distribution and also as time becomes large, the 
transient response will approach an exact 
steady-state solution. A unique feature of the 
temperature interpolation presented in equation 
(13) is that the temperature Interpolation within 
an element is time dependent. Typical behavior of 
conventional, nodeless parameter and nodeless 
variable interpolation functions during a 
transient response are compared in Fig. 2. The 
figure shows that the conventional and nodeless 
parameter temperature Interpolations retain the 
same shape during the response, but the amplitude 
of the nodeless variable temperature interpolation 
varies throughout the response. 

Rod Element 


Characteristics of the thermal and structural 
models of a rod element are shown in Fig. 3. Four 
heat transfer cases (Table 1) are considered in 
the thermal model with axial conduction combined 
with internal heat generation (source or sink), 
surface heat flux or surface convection. For 
transient heat transfer, an energy balance on a 
small segment of the rod gives the governing 
differential equations for the temperature 
distribution, T(x,t); 
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3x^ 
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pcA^ + hpT 

= hpT 

at 



(Case 1) 
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(15b) 

(Case 3) 

(15c) 

(Case 4) 

(15d) 


where A is the element cross-sectional area and p 
is the perimeter. The temperature Interpolation 
functions are first derived by solving for 
steady-state solutions to the differential 
equations above. The nodeless parameters and 
exact temperature interpolation functions are then 
written in the form of equation (10); the results 
are shown in Table 2. For Case 1 the results are 
the same as the conventional F.E.; linear 
interpolation functions yield the exact solution. 
For Cases 2-3, the results are modified by the 
Nq(x) term which introduces a parabolic variation 
of temperature from the nonhomogeneous terms, the 
heat loads. For Case 4, a completely new form for 
the interpolation function is obtained which heat 
transfer analysts recognize as the solution for a 
fin with prescribed end temperatures. 


The exact temperature interpolation functions 
(Table 2) are used to derive the element 
capacitance and conductance matrices and heat load 
vectors by evaluating the integral definitions 
given in equation (3); the results are presented 
in Appendix A. 


Use of the exact rod element interpolation 
functions (Table 2) and the element matrices yield 
exact steady-state rod temperatures, but the full 
benefit of the results is not realized unless 
Incorporated consistently in the structural 
analysis. 


For the structural response a force balance 
on a small segment of the rod (Fig. 2) gives the 
governing differential equation for the member 
axial displacement, u(x,t) 

EA 32u/8x 2 = aEA sT/ax (16) 


where E is the modulus of elasticity and a is the 
coefficient of thermal expansion. In the 
conventional two node structural element the 
displacement varies linearly and yields an exact 
stiffness matrix. For mechanical loads 
applied at the nodes, exact nodal displacements 
are obtained. Yet, in conventional 
thermal -structural analysis a linear variation of 
temperature is assumed in computing the thermal 
forces from equation (6b). Since a linear 
variation is an approximation to the true 
temperature distribution for Cases 2-4, only 
approximate thermal forces are obtained thereby 
reducing the accuracy of the structural analysis. 


In the integrated thermal -structural analysis 
the temperatures (Eq. 10 and Table 2) are employed 
to derive a consistent set of thermal forces. 
First, the temperature interpolation functions 
given in Table 2 are differentiated and 
substituted into the right-hand side of Eq. (16) 
and the differential equations are solved. The 


5 



exact displacement Interpolation functions are 
derived in the form of Eq. (10) and are shown in 
Table 3. The results for Ni(x) and N2(x) are the 
same for all four cases and are the conventional 
linear interpolation functions. The conventional 
linear tenperature distribution is modified by the 
Nq(x) term yielding in each case a nonlinear 
displacement variation within the element. The 
nonlinear variation depends on the element nodal 
temperatures and hence the thermal loading. 

The exact temperature Interpolation functions 
(Table 2) and the exact displacement Interpolation 
functions (Table 3) are used to derive the element 
stiffness matrices and equivalent nodal forces by 
evaluating the Integral definitions given in 
equation (6); the results are presented in 
Appendix A. 

The element stiffness matrix is the same as 
for conventional elements, but the force vectors 
differ from the conventional forces. These nodal 
forces, when computed using the exact nodal 
temperatures from the thermal analysis, will give 
exact nodal displacements for all heat load cases 
considered In contrast to the conventional 
approach which gives exact nodal displacements 
only for Case 1. Values for displacements within 
an element may also be computed if desired using 
Eq. (10) and the interpolation functions given In 
Table 3. 

Axl symmetric Element 

Characteristics of the thermal and structural 
models of a one dimensional axisynnetric element 
are shown in Fig. 4. The thermal model considers 
radial conduction combined with internal heat 
generation. Specified surface heating or surface 
convection on the inner and outer cylinder 
surfaces are considered through the boundary 
conditions. The governing differential 
equatlons^^ for the temperature T(r,t) are 

-k Ip (rf ) + pc = 0 (Case 1) (17a) 

-k p Ip (r|p) + pc 1^ = Q (Case 2) (17b) 

The corresponding nodeless parameters and element 
interpolation functions are derived as described 
previously and are shown In Table 4. For both 
cases the logarithmic temperature interpolation 
functions differ significantly from conventional 
linear Interpolation functions. For surface 
boundary conditions of specified heating or 
convection the same Interpolation functions 
produce exact steady-state temperature 
distributions. Element conductance and 
capacitance matrices and heat load vectors for 
Cases 1 and 2 and surface boundary conditions are 
presented in Appendix B. 

For the structural response, a 
two-dimensional plane stress or plane strain 
elasticity model is assumed (Fig.. 4). For plane 
stress the governing differential equatlon^^ fQp 
the radial displacement u(r,t) is 



where v is Poisson's ratio. The formulation of 
the plane strain problem Is mathematically 


analagous; therefore, only the equations for the 
plane stress formulation are presented herein. 

The corresponding equations for plane strain may 
be obtained by substituting equivalent elastic 
constants^^ for E, v and a. The temperature 
gradient Interpolation functions (Table 4) are 
substituted into the right-hand side of Eq. 18, 
then the differential equations are solved and 
written in the form of Eq. (10). The element 
displacement interpolation functions are shown in 
Table 5 written in terms of the nodeless parameter 
Tq. The element stiffness matrix and equivalent 
nodal forces are presented in Appendix B. 

For the steady-state case, the element 
stiffness matrix and the load vector computed 
using the exact temperature solution produce exact 
nodal displacements. Values of the displacements 
and stresses computed within the element are also 
exact. In contrast, conventional axi symmetric 
elements predict approximate element displacements 
and stresses. Moreover, such conventional 
elements typically predict stresses with lower 
accuracy than displacements. For the transient 
case, the integrated axisymmetric element predicts 
approximate displacements and stresses due to 
errors in the transient temperatures. 

Conventional elements also predict approximate 
transient displacements and stresses but with loss 
of accuracy due to temperature errors and the 
approximate formulation of the structural element. 


Thermal -Structural Analysis Programs 

The computer programs used to evaluate the 
integrated transient thermal -structural analysis 
approach are briefly described. 

TAP-STAP 

Integrated finite element thermal and finite 
element structural analysis is performed by 
linking TAP3 and STAP. TAP3 is an exploratory 
thermal analysis program similar to an earlier 
program TAP2^^ for steady-state and transient 
thermal analysis of convectively cooled 
structures. TAP3 is currently being used to 
develop thermal element methodology for Integrated 
thermal -structural analysis. STA^^ is 
an educational finite element program which was 
modified for this study to Include the Integrated 
rod and axisymmetric structural elements. TAP3 
and STAP are linked through common mass storage 
files. For the rod element, element forces per Ea 
are transmitted from TAP3 to STAP. For the 
axisymmetric element, element nodeless variables 
and nodal temperatures are transmitted from TAP3 
to STAP. 

SPAR 

SPAR^5>^^ a general purpose finite element 
program developed for production-type structural 
analysis and recently extended for thermal 
analysis was used for the conventional 
thermal -structural analyses. The SPAR program 
consists of processors which communicate through a 
data base to perform basic analysis tasks as shown 
schematically In Fig. 1(b). The incorporation of 
Integrated thermal -structural elements Into SPAR 
would yield a fully Integrated thermal -structural 
analysis capability as depicted In Fig. 1(b). 
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Appi ications 

The Integrated thermal -structural elements 
are applied to three transient examples of 
increasing complexity: (1) a rod with surface 

convection, (2) a cylinder with internal heating, 
and (3) a wind tunnel ceramic nozzle insert. The 
applications demonstrate relative performance of 
conventional and integrated finite elements. In 
addition, applications (1) and (3) compare finite 
element solutions with exact analytical 
solutions. All transient finite element 
temperature computations are made with implicit 
time-marching algorithms. 

Rod with Surface Convection 

A rod with surface convection and specified 
end temperatures is shown in Fig. 5a. Initially 
the rod is in thermal equilibrium being cooled by 
convection to a medium with T„ = 255 K, but at 
t = O'*' the convective exchange temperature is 
raised Instantaneously to T„ = 589 K thereafter 
heating the rod. A thermal transient ensues with 
the rod temperature approaching a new equilibrium 
state for large time values. Rod temperatures are 
computed from: (1) the exact analytical solution 

to equation 15d, (2) a conventional finite element 
model, (3) the nodeless parameter approach based 
on equation (11), and (4) the nodeless variable 
approach based on equation (13). In each finite 
element analysis two equal -length elements were 
used yielding an unknown at the center of the 
rod. Temperature distributions for t = 0, 0.01 
and 0.3 s are shown in Fig. (5b-5d). 

At t = 0, Fig. 5b, the nodeless parameter and 
nodeless variable approach predict the exact 
steady-state temperature distributions indicated 
by the solid line. The conventional finite 
element predicts the center nodal temperature 
well, but the nonlinear rod temperature 
distribution is only approximated by the 
linear temperature distribution of the 
conventional element. At t = 0.01 s, Fig. 5c, 
typical differences in the transient temperature 
distributions predicted by the three finite 
element approaches are demonstrated clearly. The 
conventional element yields a fair approximation 
to the center temperature, but only crudely 
approximates the true temperature distribution. 

The nodeless parameter element also yields a fair 
approximation to the center temperature, but 
predicts an extremely poor temperature 
distribution elsewhere. The nodeless variable 
approach gives the best approximation for the 
center temperature and an excellent representation 
for the rod temperature distribution. As the rod 
temperatures approach a new steady-state 
distribution at t = 0.3 s (Fig. 5d) the 
conventional element yields a fair approximation 
to the center temperature and crudely approximates 
the temperature distribution. The nodeless 
parameter and nodeless variable solutions give 
excellent approximations for the center 
temperature and the entire temperature 
distribution at this time. 

The performance of the nodeless parameter 
solutions at t = 0.01 s (Fig. 5b) is typical of 
the behavior of this approach at other 
intermediate times (not shown). The nodeless 
parameter approach employs a steady-state 
temperature distribution to approximate the 


transient behavior. Each element takes on a 
Steady-State temperature distribution as indicated 
by the humps in the predicted temperatures in 
Fig. 5b. These results are clearly not 
acceptable, and therefore this approach should not 
be used for transient response predictions. 
Instead, the nodeless variable approach should be 
employed since it gives accuracy superior to the 
conventional element throughout the response and 
predicts exact steady-state distributions. 

The use of the nodeless variable predicted 
temperatures in the corresponding structures 
problem produces more accurate displacement and 
stress distributions than conventional finite 
elements.^ For brevity, structural 
calculations with the rod element are omitted 
herein. The benefits of the integrated thermal 
structural analysis are clearly demonstrated in 
the following two examples with the axi symmetric 
element. 

Gylinder with Internal Heating 

The cross-section of an Infinite cylinder 
with axisyinmetric internal heating and equal 
specified surface temperatures is shown in Fig. 

6a. Initially the cylinder is in thermal 
equilibrium being internally heated by uniform 
internal heat generation Qq over one-half the wall 
thickness. At t = 0+, the internal heating is 
reduced to zero, and a thermal transient ensues 
with the cylinder temperatures decaying to a 
uniform distribution for large time values. The 
problem Is a rough approximation to a cylinder 
made from a semi-transparent material subject to 
Intense Internal radiation.^' 

Cylinder temperature distributions are 
computed from; (1) a conventional finite element 
thermal model employing twenty uniformly spaced 
solid elements, and (2) a nodeless variable finite 
element thermal model with two equal -length 
axi symmetric elements. Fig. 6b shows the computed 
temperatures at t = 0, 100 and 1000 s. At t = 0, 
the two nodeless variable elements predict the 
exact temperature distribution, and the comparison 
indicates the twenty element conventional model is 
sufficiently refined to represent the exact 
solution with negligible error. For the two 
subsequent times in the transient response, the 
temperatures are in excellent agreement indicating 
the capability of the nodeless variable approach 
to predict accurate temperatures with a smaller 
number of elements than the conventional approach. 

Transient displacements and stresses are 
computed from conventional and integrated 
structural models corresponding to the previous 
thermal models and a conventional finite element 
structural model with two equal -length solid 
elements. Figs. 6b-6c show the computed 
displacements and stresses, respectively, at times 
corresponding to the thermal analysis. At t = 0, 
exact displacements and stresses are predicted by 
the two integrated elements, and the validity of 
the twenty element conventional structural model, 
is verified by comparison with the exact 
solution. For the two subsequent times in the 
transient analysis, the two element nodeless 
variable approach yields excellent agreement with 
the displacements and stresses from the refined 
conventional element model demonstrating the 
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effectiveness of the Integrated thermal -structural 
analysis technique. 

The importance of integrating the thermal and 
structural analysis is demonstrated by the results 
from the two element structural model shown in 
Figs. 6b-6c. This model predicts displacements 
and stresses consistently too low and is 
inadequate to represent the structural response. 
The significant point demonstrated by the analysis 
is that even though the correct nodal temperatures 
were employed, the conventional two element model 
was inadequate to represent the structural 
response because the element thermal forces are 
computed from a linear temperature variation which 
is a poor approximation for the actual 
temperatures, and the conventional two element 
structural model was inadequate to represent the 
cylinder stiffness. Thus for improved 
thermal -structural analysis it is not sufficient 
to use nodal temperatures from a more accurate 
temperature calculation; the structural model must 
be consistently formulated to effectively utilize 
the improved temperature distribution and 
correctly model the structural stiffness. 

Wind Tunnel Nozzle Insert 

In a recent paper^^ ceramic nozzle inserts 
(Fig. 7a) for the Langley 8-Foot High-Temperature 
Structures Tunnel (8* HTST) were investigated in a 
study of concepts to improve nozzle throat life 
and tunnel performance. The Langley 8' HTST 
provides realistic temperature simulation of Mach 
7 hypersonic flight at altitudes, of 24 to 40 km. 
The products of methane-air combustion are used as 
the test medium, and the combustion products are 
accelerated through the nozzle throat to produce 
hypersonic flow in the test section. During 
tunnel start-up and shut-down the nozzle throat is 
subjected to severe temperature changes. The 
severe thermal environment causes cracking in 
metallic nozzle throats and led to the 
consideration of ceramic nozzle insert concepts. 

The nozzle throat insert problem serves as a 
convenient and demanding practical application for 
the evaluation of the integrated 
thermal -structural analysis approach. In Ref. 17, 
axisymmetric finite element thermal -structural 
models were employed, and the finite element model 
of the nozzle insert shown in Fig. 7a consisted of 
117 solid elements. The analysis showed that 
temperature gradients in the axial direction are 
relatively mild,* but the temperature gradients in 
the radial direction are severe. Herein, 
three-dimensional conventional and one-dimensional 
integrated axisymmetric element models (Fig. 7b) 
are employed to represent the radial 
thermal -structural behavior of the nozzle insert. 

A radial variation in the element size is used to 
better represent the thermal -structural response 
near the inner surface. The mesh spacing is the 
same as used in Ref. 17 and is identical for both 
conventional and integrated models. Temperature 
dependent thermal properties were considered in 
Ref. 17, but herein thermal properties are assumed 
constant. Plane strain (axial displacements equal 
to zero) is assumed for the structural analysis. 

A detailed description of five nozzle 
presssure and thermal load cases is presented in 
Ref. 17. Herein, only one thermal load case is 
considered, and a simplified representation of the 


heating and cooling from the flowing gas stream is 
employed. The nozzle heating is represented as a 
convective boundary condition with time dependent 
properties: 


T.(t) 

= T. 

0 < t 


= 0.2 T. 

125 < t 

h(t) 

= h. 

0 < t 


= 0.1 h. 

125 < t 


< 125 s 

< 150 s 

< 125 s 

< 150 s 


This representation assumes an instantaneous rise 
and fall of the heating during tunnel combustion 
start-up at t = 0 and combustion shut-down at 
t = 125s. Temperature histories computed by the 
nodeless variable approach for five radial 
locations on the nozzle insert are shown in Fig. 

8. These histories show that initially the 
temperature of the heated surface (point 1) 
rapidly approaches the convective exchange 
temperature, maintains nearly a constant value 
until 125 s and then declines rapidly as the 
nozzle is cooled by flow at the lower gas 
temperature. Temperatures in a thin layer near 
the surface tend to follow this trend also, but 
the outer surface (point 5) shows no change. 

Radial temperature variations at t = 1 s and t = 
126 s are shown in Fig. 9a. These variations show 
the very sharp radial temperature gradients at 
early and late times in the response, and that the 
variations are confined to the inner region of the 
insert. The corresponding variation of the 
circumferential (hoop) stress is shown in Fig. 

9b. The initial high temperature on the inner 
surface causes a sharply varying compressive 
stress in the inner region. As the temperature of 
the surface is reduced after 125 s, the inner 
surface stress drops significantly but the inner 
region maintains significant compressive stresses 
due to high internal temperatures. The portion of 
the nozzle outside of this region experiences a 
relatively small tensile stress throughout the 
response. 

Temperature distributions computed by 
integrated elements, conventional elements and an 
exact analytical solution^^ are compared at t = 1 
s in Fig. 10a. For the conventional finite 
element model, temperatures are computed for 
consistent and lumped capacitance matrices. The 
temperatures computed by the nodeless variable 
integrated element show the best agreement with 
the exact solution. The conventional elements 
with lumped capacitance predict temperatures that 
are too high, and the straight line interpolation 
gives a poor approximation to the exact 
distribution. The conventional elements with a 
consistent capacitance matrix predict temperatures 
for r/b = 0.575 that unrealistically oscillate 
below and above (not shown) the initial uniform 
cylinder temperature. An interpretation'^^ of this 
unrealistic behavior is that the oscillations are 
an indication of the need for mesh refinement. 
However, the results show that although the mesh 
of conventional elements is inadequate to 
accurately represent the temperature distribution 
at this time, the nodeless variable integrated 
elements give good to excellent results with the 
same mesh. 


Displacement and stress variations from the 
integrated and conventional finite element models 


8 



are compared v/lth an exact solutlon^^ jn Figs. 10b 
and 10c, respectively. Displacements computed 
from the integrated finite element approach show 
excellent agreement with the exact solution, but 
the conventional approach predicts displacements 
that are consistently too large. Stresses 
computed from the integrated elements also show 
excellent agreement with the exact stress 
solution, but the conventional elements also 
predict the stresses well, particularly the 
maximum stress which occurs on the inner surface. 
The reason that all methods show such good 
agreement is that the maximum stress at early 
response times is controlled by the temperature of 
the inner surface not the interior temperature 
distribution. In fact, for early response times, 
the maximum stress can be calculated^^ with good 
accuracy by 

where is the initial uniform temperature. 

Temperature and stress distributions computed 
by integrated and conventional elements are 
compared at t = 126 s in Figs. 11a and 11b, 
respectively. The figures show that the mesh is 
not sufficiently refined for the conventional 
elements to accurately represent the large 
variation of temperature and stress which occurs 
near the inner surface. The highly nonlinear 
variation of temperature and stress in this region 
cannot be realistically represented by the 
conventional element linear interpolation 
functions. The sharp discontinuities of the 
slopes of the temperature and stress curves at r/b 
= 0.575 are an indication of the inadequacy of the 
conventional element mesh. The integrated 
elements, however, predict realistic variation of 
the temperatures and stresses with good continuity 
of slopes of the temperature and stress curves. 

In comparison to the integrated analysis, the 
conventional approach with lumped capacitance 
underestimates the maximum temperature by 5 
percent and underestimates the maximum stress by 9 
percent and predicts the wrong locations for these 
maximums. 

The examples demonstrate the potential of the 
integrated approach for the one dimensional 
examples considered. However, evaluation of the 
full potential of the approach requires additional 
investigation by; (1) applications to more 
complex structures with other finite elements, (2) 
consideration of nonlinear problems due to 
temperature dependent material properties, and (3) 
consideration of other heat transfer modes 
including radiation. 


Concluding Remarks 

An integrated transient thermal -structural 
finite element analysis approach is presented. 

The integrated thermal -structural analysis is 
characterized by: (1) thermal and structural 

finite elements formulated with a common geometric 
discretization with each element formulated to 
suit the needs of its respective analysis, (2) 
thermal and structural finite elements fully 
compatible during the coupling of the two 
analyses, and (3) use of equivalent thermal forces 
based on a consistent finite element formulation. 


Integrated thermal -structural rod and one 
dimensional axisymmetric elements are developed 
and used to demonstrate the integrated analysis 
approach for conduction with surface convection. 

A unique feature of the elements is the use of 
nodeless variable interpolation functions for 
element temperatures and displacements. For 
steady-state linear analysis, the elements produce 
exact values of temperatures, displacements and 
stresses. This feature means the transient 
response can start from an exact initial 
temperature distribution, or if the response 
approaches thermal equilibrium for large time, 
exact steady-state temperatures, displacements and 
stresses are predicted. 

The integrated thermal structural elements 
are evaluated by solving three transient examples 
of increasing complexity: (1) a rod with surface 
convection, (2) a cylinder with internal heating, 
and (3) a wind tunnel ceramic nozzle insert. 
Temperatures, displacements and stresses predicted 
by the integrated approach are compared with 
results from conventional finite elements and 
exact analytical solutions. In the examples, the 
integrated approach showed a clear superiority to 
the conventional finite element approach. For the 
same number of elements, the integrated approach 
has superior accuracy. For equivalent accuracy, 
fewer integrated elements are required; for 
instance, in the analysis of the cylinder with 
internal heating, two integrated elements gave 
results equivalent to predictions from twenty 
conventional elements. The nonlinear variation of 
temperatures, displacements and stresses within 
the integrated elements permits an accurate 
representation of these distributions with only a 
few elements. This feature was clearly 
demonstrated in the nozzle insert example which 
experienced sharp radial variations of the 
temperature and stresses due to instantaneous 
application and removal of convective heating. 

The examples demonstrate the approach offers 
significant potential and should be investigated 
further by applications of other element types to 
more complex thermal -structural problems. 
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Appendix A: Rod Element 

The rod element matrices are derived from the 
general element equations (3) and (6) using the 
element temperature and displacement 
Interpolation functions given in Tables 2 and 3, 
respectively. 


Capacitance Matrices 

Cqo = pcAL/30 

^01 “ ^02 " pcAL/12 

Cii = C 22 - pcAL/3 

Ci2 ~ ^21 " pcAL/6 


(Cases 1-3) 


Cqo 

Cqi =* ^02 

Cii * C22 
C12 = C21 


pcA 

pcA 

pcA 

pcA 


r /cosh mL-l V L + J 

^ sinh mL y\sinh mL m/ J 

/ /^ 

[ sinh mL cosh mL-mL 1 
2m sinh^ mL J 


2m sinh mL 

sinh mL cosh mL-mL 
— 2 

2m sinh mL 

‘ mL cosh mL-sinh mL 
2 

2m sinh mL 


(Case 4) 


where 

m = ^hp/kA* 
Conductance Matrices 
Kqo = kA/3L 

Kqi = Kq2 = 0 

Kii = <22 = kA/L 

<12 = 1^21 = -kA/L 


(Cases 1-3) 


KOI 

Kll 

<12 


Kqo 



2(cosh mL -1) 
sinh mL 


•^02 

K22 

K21 


0 

m 

h£ 

m 


cosh mL 
sinh mL 

_ 1 

sinh mL 


•] 


(Case 4) 


For case 4, the conductance matrix Includes the 
contributions from conduction and convection given 
in equations (3b-3c). 


Heat Load Vectors 

{qF = QAL [1/6 1/2 1/2] (Case 2) 
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JqT = tlPl- [1/6 1/2 1/2] (Case 3} <^22 = (1 - 2w ♦ 2v^) - a^J 

% ■''p^[[-|(^TT?rr)] 

Ci-Cc ■>'p/.[i(^TraTf^)] 


(Case 4) 


Stiffness Matrix 

Kii = K 22 = AE/L 
Ki2 = K 21 = -AE/L 
Force Vector 


where 


(Cases 1-4) 


where w = 1 n(b/a). 

Conductance Matrix 

Kqo = kw[w(l-(a/b)4) - (l-(a/b)2)2j 

•^01 = Kq2 = 0 

Kii = K22 = k/w 
Ki 2 = ^ 2 ]^ = -k/w 

For surface convection on the cylinder inner and 
outer surfaces additional conductances are added 


and 


{F|^ = C-F F] 


to the system conductance matrix: 


e 



r = a, = ha 


F = aEACT^/S + (Ji + T2)/2] 

(Cases 1-3) 


r = b, = hb 




Heat Load 

Vectors 


F = oEACC^Tq + C2(T]^ + T 2 )] 

(Case 4) 

Qo 


■?’ J 

P _ , 2 (cosh mL -1) 

" mL sinh mL 

P _ cosh mL - 1 
^2 mL sinh mL 


Qi = 
Q2 = 

w [■ I" " * T 

(Case 2) 


Appendix B: Ax i symmetric Element 

The axi symmetric element matrices are derived 
from the general element equation (3) and (6) 
using the element temperature and displacement 
interpolation functions given in Tables 4 and 5, 
respectively. Matrices are the same for all cases 
except where noted. 


Capacitance Matrix 

Coo = pc { |^4w2(aWbW) + 9w(a'^-b'*' 


For the surface convection and surface heating on 
the cylinder inner and outer surfaces additional 
nodal heat loads are added to the system heat load 
vector: 

r = a, Q = haT^ + qa 

r = b, Q = hbT„ + qb 

Stiffness Matrix (Plane stress) 

(b^ + a^) - v(b^ - a^)J 


24 b^ 


. 6(aV'2 




Coi " - Ua^v<^-hv(7a^+3b^)(a^-b^) 

•16b w L 

+ 4(a^-b^)^J 

C02 ■ ^^9 r 4b^w^ - w 

16b‘^w L 


Kll = 

E 

1 

(1-v^) 

(b^-a^) 

Ki2 = 

E 

1 

(l-v^) 

(b^-a2) 

K 22 = 

E 

1 

(1-v^) 

(b2-a2) 


(b^ + a^) + v(b^ - a^)J 


(7b^+3a^)(b^-a^) 


* 4(aV)2] 


Force Vector (Plane stress) 

TT-!-! 

(1-v) w(b -a^) i L 

+ (2b2w-b2+a2)(T2+wT^) - w2 T, 


2a^w+b^-a^)(Tj+^ wT^) 


Cii = 2w + 2w^)j 

£C_ r 2 2 2 2 ] 

C 12 “ C 21 “ 2 M ^ w(a + b ) J 

4w^ L -I 


. _ Ea 
F2 


- 2(b2-a2)w T^j | 
^ j ^ r(-2aVb2-a2)(Ti+^ wTq) 

-2.1 2 L 9 


(1-v) w(b -a )' “• b' 

+ (2b2w-b2+a2){Tp+wT ) - w2 T 


- 2(b2-a2)w Tp j 
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Table 1 Rod Element Heat Transfer Cases 


CASE HEAT TRANSFER MODE 

1 CONDUCTION WITH NO HEATING 

2 CONDUCTION WITH INTERNAL HEAT GENERATION 

3 CONDUCTION WITH SURFACE HEAT FLUX 

4 CONDUCTION WITH SURFACE CONVECTION 


Table 2 Rod Element Interpolation Functions 


Case 

Nodeless 
Parameter, Tq 

No(x) 

Ni(x) 

N2(x) 


1 

0 

r (1- f) 

1-u 

X 

L 


2 


r (1- f) 


X 

L 


3 

a£ti 

2kA 

r (1- r) 

>-r 

X 

L 


4 

T 

CD 

- sinh m(L-x) sinh mx 
sinh ml " sinh ml 

sinh m(L-x) 
sinh ril 

sinh mx 
sinh ml 



WHERE m 


Table 3 Rod Element Displacement Interpolation Functions, N 5 (X)* 
CASE Nq (X) 


1 


2 


(x2 


X) 


2 

3 


a(T2-Tj) L 
2 

a(T2-Tj) L 
2 


(X - X) + — (- X + 3X"^ - 2X'^) 

n qT L « o 

- X) + — ^ (- X + 3X"^ - 2X*^) 


' ^2“^l ^ (cosh ml -1) r 1 

4 a m sT nh lil ^ 

*For all cases. N^IX) = 1 - X N2(^F= X“" 
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Table 4 Axisymmetric Element Temperature 
Interpolation Functions 


Case Nodeless Parameter, Tq 

1 0 
2 Qb2/(4kw) 

N(,(r) = 1n(r/a) + a2in(b/r)/b2 . r^w/b^ 

Ni(r) = ln(b/r)/w 
N 2 (r) = ln(r/a)/w 

where w = ln(b/a) 


Table 5 Axisymmetric Element Displacement 
Interpolation Functions 





X 


(a) Conventional. 



X 


(b) Nodeless parameter. 



X 


(c) Nodeless variable. 

Fig. 2 One dimensional element interpolation 
functions. 



• THERMAL AND STRUCTURAL • IMPROVED THERMAL ELEMENTS 

ELEMENTS SIMILAR - FORMULATION FUNCTION OF HEATING 

• NODAL TEMPERATURES jTl • COMPATIBLE THERMAL DATA TRANSFER 

TRANSFERRED AS REQUIRED BY STRUCTURAL ELEMENT 

• THERMAL FORCES BASED • THERMAL FORCES BASED ON ACTUAL 

ONLY ON NODAL |T[ TEMPERATURE DISTRIBUTIONS 


(a) Conventional analysis, (b) Integrated analysis. 


Fig. 1 Conventional versus integrated thermal and 
structural analysis. 



GENERATION HEAT FLUX 
Q (t) q It) 

THERMAL MODEL 

Fig. 3 Thermal and stress models of rod element. 
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ONE DIMENSIONAL THERMAL MODEL 


Fig. 4 Thermal and stress models of ax 1 symmetric 
element. 
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(b) Comparative temperature distributions at t = Os. 
Fig. 5 Continued. 


255 K 


h.Too 


T=31l K 



T = 533 K 


(a) Rod heated by surface convection. 


Fig. 5 Conventional and integrated finite element 
solutions for a rod with surface 
convection. 
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(c) Comparative temperature distributions at t = .01s. 


Fig. 5 Continued. 
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x/L 

(d) Comparative temperature distributions at t = .30s. 
Fig. 5 Concluded. 



(a) Cylinder with internal heat generation. 

Fig. 6 Conventional and integrated finite element 
solutions for cylinder with Internal 
heating. 
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(b) Comparative temperature distributions. 
Fig. 6 Continued. 
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(c) Comparative displacement distributions. 
Fig. 6 Continued. 
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(d) Comparative circumferential stress distributions. 
Fig. 6 Concluded. 
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(a) Nozzle temperature distributions. 

Fig. 9 Nozzle Insert temperature and stress 
distribution computed by integrated 
thermal -structural finite elements. 


(a) Nozzle cross section. 



MODEL MODEL 

(b) Finite element models at section A- A. 

Fig. 7 Tunnel nozzle Insert and 

thermal -structural models. 



Fig. 9 Concluded. 
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(a) Comparative temperature distributions. 

Fig. 10 Comparison of nozzle temperatures, 

stresses and displacements at t = Is. 
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(b) Comparative displacement distributions. 
Fig. 10 Continued. 
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(c) Comparative circumferential stress distributions 
Fig. 10 Concluded. 
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(a) Comparative temperature distributions. 

Fig. 11 Comparison of nozzle temperatures and 
stresses at t = 126s. 



.56 .57 .58 .59 


r/b 

(b) Comparative circumferential stress distributions. 
Fig. 11 Concluded. 
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